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Abstract 

The large rate of the J/ijj-\-rj c production observed by the Belle collaboration has posed a serious 
challenge to our understanding of the quarkonium. We examine a scenario that there exists a light 
scalar boson around 5 GeV, with an enhanced coupling solely to the charm quark. Such a scenario 
would explain the Belle anomaly. It also predicts a significant increase in J/ip plus open charm 
pair production. An immediate test for the scalar boson is to look for a peak around 5 GeV in the 
recoil mass spectrum of J/ip, because the scalar boson can also be produced in association with 
the J/ip. Finally, we also point out that the process e + e~ — ► H'y is sizable for observation. 
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I. INTRODUCTION 



The best understood QCD system is perhaps the heavy quarkonium, made up of heavy 
quark and antiquark. Known systems include charmonium, bottomonium, and B c mesons. 
We thought we understand these systems, but surprises often come up as more and more data 
become available. A well-known incident was the surplus of J/ip produced at high transverse 



momentum at the Tevatron 111]. The color-singlet model 31 missed the cross section by 

n 

more than an order of magnitude. The color-octet mechanism [3] and subsequently the 



nonrelativistic QCD model |4| were formulated to explain the surplus. 

Recently, the Belle Collaboration observed double charmonium production, namely, 
J/ip + rj c at y/s = 10.6 GeV |5|. The production rate is anomalous 

sag 

<r(e + e~ — > J/ip rjc) x B(j] c — > 2 charged) 

= (0-033 t° Z ±0.009) pb, 



which is roughly an order of magnitude larger than any QCD predictions [8|]. It poses 
a strong threat to our understanding of quarkonium. However, one has to note that the 
J/ip is detected via its leptonic decay while r] c is inferred from the recoil mass spectrum 
of the J/ip. Therefore, the r] c is actually not identified, and perhaps it could be a J/ip, or 
something else. Based on this argument, Bodwin, Lee, and Braaten [9J argued that Belle 
in fact observed double J/ip*s, which are produced via double virtual photons. However, 
the further analysis by Belle does not seem to support this idea p. Thus, the anomalous 
production remains one of the challenging problems in QCD, and arouses speculations. For 
example, a heavy glueball of a mass very near m Vc with a rather large matrix amplitude has 



been proposed [lOj to account for the excess. 

In this Letter, we study a possible scenario of the existence of a light scalar boson H, 
which has an enhanced coupling to the charm quarks. It contributes to both the J/ip + r\ c 
production, and J/ip plus open charm pair production. We show the allowed parameter 
space region of the mass of this scalar and the coupling to charm, fitted to the Belle data. 
We find that a coupling constant gjj w 1 and tuh ~ 5 GeV can fit the data well. An 
immediate test for this light scalar boson is to look for a peak around 5 GeV in the recoil 
mass spectrum of J/ip, because this scalar can also be produced directly in association with 
a J/ip. So far, the Belle collaboration has only shown the recoil mass spectrum with very 
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FIG. 1: One of the four contributing Feynman diagrams for e + e — > J/iprj c via an intermediate 
scalar boson H. 

fine bin size up to 4 GeV sag , where ?7 C , Xco, and r) c (2S) can be seen. However, the bin 
size in the region above 4 GeV is too large. Therefore, we urge the Belle collaboration to 
scan with very fine steps in the region above 4 GeV. (Even with a large bin size, we notice 
that there may be something peculiar at around 5.2-5.4 GeV bins j^, |fj.) Finally, we point 
out that another process e + e~ — > H^f, which is negligible within the SM, becomes sizable in 
the our scenario. This is another interesting test for the scenario. 



II. e + e — ► J/ipri c 

Suppose the interaction of the scalar boson H with the charm-quark pair is given by 

C = -g H ccH, (1) 

with a coupling constant gn that will be determined later. 

Assuming that the scalar decays into a cc pair only, the width of H is given by 

3/2 



3 2 



4m c 

m H 



(2) 



For m,u ~ 5 GeV and Qh ~ 1, the r# ~ 0.06m#, which is narrow compared with its mass 
but much wider than the widths of J/ij) and r\ c . In the following, we use the Breit-Wigner 
prescription for the scalar boson propagator. 

One of the contributing Feynman diagrams for e + e _ — > J/ipr] c is shown in Fig. CI 
The sum of the amplitudes of the four contributing Feynman diagrams for e~(pi)e + (p2) — > 
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2„2^ ^ ^(0)^(0) 1 



J/i[)(ki)r] c (k2) is given by 

iM = -4g n e'Q e Q c , , , 

7r s z s/4 — m H + il h^tlh 

x v(p 2 )Yu(p 1 )e fipXa k p 2 k^e\k 1 ) (3) 

where e(fci) is the polarization 4- vector of the J/ip, and s is the square of the center-of- 
mass energy of the collision. it^(O) and -Rr/ c (0) are the wavef unctions at the origin, with 
the normalization J °° \R{r)\ 2 r 2 dr = 1. Here the center-of-mass energy ^/i of the collision 
is at 10.6 GeV, and so we simply only take the intermediate photon approximation (the 
intermediate Z contribution is suppressed by s/m 2 z ~ 0.013.) 
The angular distribution of the process is given by 



dcos9 e c it s 2 

x— — — f3 3 {I + cos 2 6) , (4) 



[s/4 - m 2 H ) 2 + T 2 H m 2 H 
where 



ntj, m 2 \ m 2 ,, m 2 




/3 = X I[1--*-^) -4^^ 
s s J s s 

Integrating the angle we obtain the total cross section 

3 7T S Z 

1 



X (s/A-m 2 H ) 2 + T 2 H m 2 H P3 (5) 

Let us examine the allowed parameter space of mu and gu that can fit the Belle cross 
section. Combining the systematic and statistical errors of the Belle data the la range is 
0.033 ± 0.011 pb. We show the contours of the central value, ±lcr, and ±2cr in Fig. |21 
Taking gg = 1 the allowed central value (2a range) for mu is 

m H = 5.10 (4.87 - 5.20) GeV or 5.45 (5.36 - 5.66) GeV . (6) 

Here we used \R^(0)\ 2 = \R V (0)\ 2 = 0.8 GeV 3 . 

It is straightforward to adapt our scenario to explain other peaks in the recoil mass 
spectrum between 3 — 3.8 GeV, i.e, Xco and rj c (2S) Q]. Since our scenario only re- 
quires the van at around y/s/2, it would also explain the enhancement to the produc- 
tion of J/i/j + (rj c (lS),Xco, r lc(2S)). Note that our scenario would not give enhancement 
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FIG. 2: The ±1<7 and ±2<r bound for the (mH,gH) plane due to the Belle data of 0.033 ± 0.011 
pb. 

to J/ip + (J/ip, Xci, Xc2, ip(2S)) because of the spin-parity-charge-conjugation of the photon 
propagator. This is in agreement with the observation by Belle 

Note that if we choose a smaller gjj, the corresponding central value m# of the fit to the 
Belle data would become closer to \fsj2 = 5.3 GeV, and the 2a range would also become 
smaller, as shown in Fig. |21 The smaller the value of gn, the more fine-tuned value for win 
that the fit gives. So in the following predictions of our scenario we simply take gu = 1 and 
the results just scale as some powers of gu- 

If we replace the scalar H by a pseudoscalar boson A , the diagram similar to Fig. 1 with 
H replaced by A gives vanishing amplitude in the static approximation. There is another 
diagram with the virtual A bremsstrahlung off the charm quark turning into the r] c state. 
However, such a contribution is important only when the mass m^o is around 3 — 4 GeV, 
which is unfavorable based on the constraint of the Belle search. 

The static approximation of using the wavefunction at the origin may not be so valid when 
the Higgs pole is involved, however our simple calculation serves the purpose of estimating 
of the size of the production rate due to the light Higgs boson H. 

III. e + e~^J/ipcc 

This process is very similar to the original J/ip + 7] c production, except that the cc pair 
does not form a bound state but goes into an open pair. They will hadronize into D mesons. 
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FIG. 3: The production cross section for e + e — > J /tpcc Vs ran- The value of gu = 1 is chosen. 

The calculation is straightforward, but however the result is not simple enough to be put 
here. We show the cross section verse m# at \Js = 10.6 GeV in Fig. El The cross section 
scales as gjj. At m H ~ 5 GeV, a(e + e~ — > J/ipcc) ~ 0.15 pb for g H = 1. Therefore, it 
contributes substantially to J/ip+ open charm production. 



IV. e+e" -» ipH 



The sum of the two Feynman diagrams of the process e G»i)e + (j?2) 
given by 

%M = 2V3e 2 Q e Q c9H i i— — 

y/Trfn^ s s -m^ + mjj 

xv(p 2 )^u{ Pl ) [(a - m 2 H + m%Y{k x ) - 2k 2 ■ e(ki)k?] 



ii(h)H(k 2 ) is 



The differential cross section is 



i 



dcosO 



777 I 

+ V(l + cos 2 
s 2 



where 



P 



ml m 



H 



-4 



s s 



The integrated cross section is 
a = Qa 2 g 2 H Q 2 e Q 



2 ,2 / -,2/ -,2 1-^(0) | 2 



(s - m| + m|-) 2 







(7) 



(8) 
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FIG. 4: The production cross section for e + e — ► J /ipH Vs m#. The value of gn = 1 is chosen. 

In Fig. El we show the production cross section of a(ipH) with gjj = 1. The cross section 
scales as g 2 H . The scalar boson H so produced is heavy enough to decay into an open charm- 
quark pair. Thus, it will give rise to open charm pair production, with the pair around 5 
GeV. If the Belle could scan the recoil mass spectrum of J / ip further up to 5 — 6 GeV, they 
should be able to see if this scalar boson H exists. This is the most immediate test for the 
scenario. 

V. e+e- -» #7 

In the SM, the leading contribution to the process e + e~ — > H'j comes from the one-loop 
diagrams with the W boson and heavy fermions running in the loop. In our scenario with H 
coupling strongly to charm, this channel could be sizable for detection. The SM calculation 
was completed in Ref. We use its formulas with the charm-quark contribution only, 

and modify the SM coupling to our enhanced coupling as follows: 

V2G F m 2 c -> g% . 

We show the ratio of cx(e + e~ — > Hj) / 'er(e + e~ — > at y/s = 10.6 GeV for m H = 3 — 9 

GeV with gjj = 1 in Fig. |SJ The ratio scales as g\. Note that with the SM coupling strength 
(gm c /2mw), this ratio is very small, of order O(10~ 9 ). With the enhanced charm coupling, 
the ratio is of order 5 x 10~ 5 for tuh = 5 GeV. Thus, the cross section of e + e~ — > H'j for 
uih = 5 GeV at y/s = 10.6 GeV is about 40 fb, which will give thousands of events at Belle 
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FIG. 5: The ratio of a(e+e' Hj)/a(e + e~ -> at yfs = 10.6 GeV for m H = 3 - 9 GeV. 

The scenario is gu = 1. 

or BaBar with an integrated luminosity of 100 fb _1 . 
VI. DISCUSSION 

One may concern about the high energy analog of e + e~ — > ccH production. In high 
energy collisions, we would expect extra cccc in the final state, with one of the cc pairs 
clustering at around 5 GeV. To our knowledge there have not been any measurement on 
four-charm production, nor any searches explicitly looking for anomalous four-charm events 
at LEP or at the Tevatron. There were perhaps some mysterious 4-jet events recorded by 
ALEPH [12], but it had no explicit identification of charm in those events. Also, the search 
for four-charm events suffers severely from the weak charm identification, which is, to some 
extent, overlapped with the b-quark identification. Therefore, there is no explicit constraints 
on four-charm production coming from high energy experiments. 

Another possible source of constraints comes from Upsilon T(IS') decays. From the 
Particle Data Book Q, the branching ratio B(T{1S) -> J /if) + X) = 1.1 ± 0.4%. Given 
such a large error it is impossible to rule out our scenario, because J /if) + r] c , J/ij) + H, or 
J/ip+ open charm could only give a branching ratio below 0.4%. There is also a limit on 
B(T — > 7 + X) < 3 x 10 -5 , where X is a pseudoscalar boson. In our scenario, we have a 
scalar boson and so the limit is not directly applicable. Even so, B(T — > Hj)/B(T — > fi + fi~) 
would be very similar to cr(e + e~ — > H'~f)/a(e + e~ — ► ~ 5 x 10~ 5 (from Fig. EJ) without 
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an explicit calculation. It implies B(T — > 7 + H) m 10~ 6 , which is well below the existing 
limit that we stated above. 

The anomalous double charmonium J /if) + r\ c production measured by Belle posed a 
very strong challenge to our understanding of QCD. We have investigated the scenario 
of a light scalar boson of mass about 5 GeV with a large coupling (of order the weak 
coupling or even larger) to the charm quark. We found that it could explain the J /if) + 
1] C anomaly, and provided the much-needed J/ip plus open charm production. It would 
also explain the enhancement to the production of J /if) + (r) c (lS), Xco, Vc(^S)), but not to 
J /if) + (J/if>, Xd, Xc2, if)(2S)) because of the spin-parity-charge-conjugation of the photon 
propagator. This is in agreement with the observation by Belle. An immediate test for this 
light scalar boson is to scan the recoil mass spectrum of J / if) and look for a peak around 5 
GeV, because this scalar can also be produced directly in association with a J /if). So far, 
the Belle collaboration has only shown the recoil mass spectrum with very fine bin size up 
to 4 GeV E3], where r] c , Xco, and r] c (2S) can be seen. However, the bin size in the region 
above 4 GeV is too large. We, therefore, urge the Belle collaboration to scan the recoil mass 
spectrum with very fine bin size above 4 GeV. (Even with a large bin size, we notice that 
there may be something peculiar at around 5.2-5.4 GeV bins 5], E|, but this is only our wild 
guess.) Finally, we have pointed out that another process e + e~ — > Hj, which is negligible 
within the SM, becomes sizable in the new scenario of a light scalar boson, which couples 
strongly to charm. This is another test for the scenario. 

One may ask where this scalar boson comes from. Such a light Higgs may come from 
the general form of the multi-Higgs doublet sector. The required enhanced coupling to the 
charm quark implies dedicated fine tuning. It also requires additional fine tuning to kill the 
unwanted flavor changing neutral processes. If we stay within the scope of the two-Higgs 
doublet sectors with the quark couplings like those in the minimal supersymmetric standard 
model, the coupling strength of H to charm is gm c cot /3/2mw- Therefore, in order to have 
an effective coupling g# ~ g we need a cot (3 ~ 2myi//m c ~ 0(100). Such a large cot (3 is 
unnatural and perhaps causes a big problem to the top Yukawa coupling. We do not offer 
a good answer to the structure of the Higgs sector in this paper. 
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